Superexchange induced canted ferromagnetism in dilute magnets 
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We argue, in contrast to recent studies, that the antiferromagnetic superexchange coupling be- 
tween nearest neighbour spins does not fully destroy the ferromagnetism in dilute magnets with 
long-ranged ferromagnetic couplings. Above a critical coupling, we find a canted ferromagnetic 
phase with unsaturated moment. We have calculated the transition temperature using a simplified 
local Random Phase Approximation procedure which accounts for the canting. For the dilute mag- 
netic semiconductors, such as GaMnAs, using ab-initio couplings allows us to predict the existence 
of a canted phase and provide an explanation to the apparent contradictions observed in experi- 
mental measurements. Finally, we have compared with previous studies that used RKKY couplings 
and reported non-ferromagnetic state when the superexchange is too strong. Even in this case the 
ferromagnetism should remain essentially stable in the form of a canted phase. 
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I. INTRODUCTION 

The physics of disordered/dilute magnetic systems has 
attracted a considerable interest and attention from both 
theoreticians and experimentalists. Among these ma- 
terials one finds for instance manganites (LaSrMnOs, 
LaCaMn03,..)fii^'^ diluted magnetic semi-conductors 
as GaMnAs^ which were widely studied, the so-called 
dP materials (Hf02, Ga,0,...)r^ the Heusler alloys as 
Ni2MnSn^i^ or the double perovkites as Sr2FeMo06'^'ii 
In these materials one of the key issue is the under- 
standing of the influence of the carrier (hole/electron) 
concentration on both magnetic and transport proper- 
ties. Indeed, the variation of the carrier concentration 
often leads to drastic changes and gives rise to interesting 
physics. In particular a competition arises between di- 
rect or superexchange™ interaction of the localized mag- 
netic moments and indirect couplings via the itinerant 
carriers. In general, the superexchange coupling dom- 
inates at low carrier concentration but is overtaken by 
the ferromagnetic contribution at higher concentration. 
For example, in manganites (non dilute) the superex- 
change coupling competes with the double exchange cou- 
pling and leads to canted ferromagnetic phasesj ^'^d'^d^ 
However, as soon as disorder is introduced into the sys- 
tem, new magnetic phases may appear such as ferromag- 
netic droplets in a canted antiferromagnetic matrix as ob- 
served in manganites li^iii In dilute systems, where the 
probability to have nearest neighbour pairs is small it 
is not clear whether the superexchange coupling has the 
same effects. In particular, it is not obvious that su- 
perexchange alone can eventually completely destroy the 
ferromagnetic phase or induce new phases. The aim of 
the present study is to focus on this issue. 

In this paper we show that in a dilute system of 
classical spins the superexchange competes with the 
long-ranged ferromagnetic couplings and favors a canted 
ferromagnetic phase in part of the phase diagram 




FIG. 1: (Color online.) Schematic representation of the 
canted ground-state resulting from the competition between 
the long range ferromagnetic couplings ,Jij and the superex- 
change coupling Jaf. The spins involved in pairs get canted 
and the angles 6i vary from spin to spin. 



(temperature-concentration). However, in contrast to 
non dilute materials and double exchange systems, only 
spins involved in nearest neighbour pairs get canted 
(Fig. [1]). This is particularly relevant for diluted mag- 
netic semiconductors (DMS). In DMS, the magnetic cou- 
plings are extended and the superexchange dominates at 
sufficiently low carrier concentration. In these materi- 
als, a conflict between the measured low temperature to- 
tal magnetic moment obtained by SQUID measurements 
and the density of spins extracted from X-ray diffrac- 
tion (XRD) is often observedJ^i^i2ii22 We shall see 
that the existence of a canted phase provides a natu- 
ral explanation to the observed disagreement. We also 
solve a conflict between recent Monte Carlo simulations — 
which found ferromagnetism in a region where the Self- 
Consistent-Local RPA predicted an instabilityi^ This in- 
stability actually signals a new phase with unsaturated 
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ferromagnetism, as we shall see. Because the nature of 
the ground state was not analyzed, this conclusion was 
missed in the Monte Carlo studies which focused on the 
amplitude of the Curie temperature only. 

The manuscript is organized as follows. In the first 
part, we will analyze the effect of the superexchange cou- 
pling assuming a simple model for the extended exchange 
integrals. In the second part we will discuss the specific 
case of GaMnAs where it is known that superexchange 
coupling dominates over the indirect ferromagnetic con- 
tribution for sufficiently low hole density. In this part, 
for a quantitative study, realistic couplings will be taken 
from ab-initio calculations (TB-LMTO).— In the third 
and last part, we will discuss the case where the cou- 
plings are of the RKKY form, in order to study the com- 
petition between superexchange and frustration effects 
induced by the oscillating tail. It will be shown that, in 
the presence of the superexchange coupling, the stabil- 
ity region is significantly larger than found in previous 
studies.^ 

In the following, we will consider the diluted Heisen- 
berg Hamiltonian which reads, 

H = — ^ ' ^ ^ XjXjJjjSj ■ Sj H- ^ ^ XjXj J AF Sj ■ Sj (1) 

where the random variable Xi is 1 if the site is occupied 
by a magnetic impurity (otherwise 0) . The total concen- 
tration of magnetic impurities is x. The localized spin 
Si at site i is classical {\Si\ — 1). The first term cor- 
responds to the long-range exchange couplings and the 
second term is the nearest neighbour antiferromagnetic 
superexchange contribution. Because we will discuss the 
particular case of GaMnAs, for convenience we have per- 
formed all the calculations for a fee lattice. 



II. A SIMPLE UNFRUSTRATED MODEL 

In this section we consider a simple model where the 
couplings are relatively extended but all ferromagnetic 
Jij = jQe~~^. The parameter A controls the range of 
the couplings; but since there is no abrupt cut-off, there 
is no strict percolation threshold for the ferromagnetism 
in this problem. The tail of the couplings always induce a 
finite transition temperature. In fact, this model is not so 
far from the exchange couplings in III-V diluted magnetic 
semiconductors such as Gai_a:Mna:As or Gai_a;Mna;N as 
calculated from first principles. More realistic couplings 
will be used in the next section. Whilst for Jaf = and 
at T = K, the ground state is ferromagnetic and fully 
saturated (no frustration), we discuss its nature in the 
presence of Jaf- For this purpose, for a given configu- 
ration of disorder (position of the magnetic impurities), 
we minimize numerically the total energy associated to 
the Hamiltonian ([1]) with respect to the angles (6*^,0;). 
For simplicity, we consider only the case where the spins 
are coplanar (0^ = 0). The calculations were performed 



for systems containing typically 1000 diluted spins. We 
have found that beyond a critical value of Jaf, one pair 
of nearest neighbour spins {Si,Sj) starts to get canted 
(see Fig. [Ij whereas all the other spins remain aligned 
along the magnetization axis. As we increase Jaf more 
pairs get canted but the spins that have no nearest neigh- 
bour remains almost parallel to the z axis. The canting 
results from the competition between the local field re- 
sulting from the long range couplings and the superex- 
change contribution of the nearest neighbour spin. Each 
spin sees a different environment and therefore the cant- 
ing does not occur simultaneously for all the pairs as we 
increase Jaf- For the same reason, the canting angles 
are different from spin to spin. For the unpaired spins, 
however, the canting angle is very small. This results 
from two combined effects. On one hand, for each pair, 
9i is close to —0j, so that their resulting transverse field 
is small. On the other hand, a given unpaired spin Sk 
experiences the sum of the transverse fields due to all 
canted pairs. Because of their random sign, the sum av- 
erages out to a small value. We can therefore neglect the 
small canting angles of the unpaired spins, as shown in 

Fig.m 

For simple illustration, we recall what happens to a 
single pair in the effective field of the other spins. The 
energy of this pair is E = Jaf cos{di — Oj) — hi cosOi — 
hjCosOj, where hi = "^iJn is the local field on spin i 
(for simplicity, we assume hi ~ hj ~ h to be the same 
for both sites, so that Oi = —6j = 6). The minimization 
gives a canting angle. 



for Jaf > h/2 (and = Q otherwise). For Jaf ~^ oo, 
the two spins are anti-aligned and orthogonal to the other 
spins: they are effectively decoupled from them. We em- 
phasize that, in our calculations, we have kept the real 
local fields hi that differ from site to site. In addition 
we also have a finite probability to have trimers of spins, 
quadrimers, etc. so that the real canting angles are not 
given by ^ but are determined self-consistently. 

Now, in order to calculate the critical temperature, 
we include thermal fluctuations about the ground state. 
In the case where the ground state is fully polarized, 
it has been shown that the self-consistent local random 
phase approximation method (SC-LRPA)^i^i2^ is reli- 
able. When directly applied to models in presence of 
superexchange couplings, it has been argued that ferro- 
magnetism disappears when the nearest neighbour cou- 
pling is dominated by the antiferromagnetic superex- 
change contribution. ^"^'^^ In fact, the reported instability 
turns out to occur at the same critical value at which the 
ground-state becomes canted, as found above. The in- 
stability therefore directly reflects the change of ground 
state and should not be interpreted as the result of frus- 
tration of the long-range couplings. In order to correct 
for this, we extend the SC-LRPA to calculate the criti- 
cal temperature of a canted state. In principle, as in the 
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original SC-LRPA one has to start with the equation of 
motion of the retarded Green's function, Gf^{uj), where, 



dte^-'9{tm\t),SU0)]) (3) 



where /x, f are the spin components. Because of the 
canted ground-state, the decouphng of the equation of 
motion involves both the longitudinal (S^) and trans- 
verse magnetizations (Sf). As a consequence, the trans- 
verse Green's function, G^~, is now coupled to both 
G^j^ and G^^ . Since solving these coupled equations 
is more involvedfSi, we propose a simplified ansatz: al- 
ter the determination of the canting angles {9i} (for a 
given configuration of disorder) we map the canted prob- 
lem to an effective fully ferromagnetic one with reduced 
spin amplitude. Si — > S cos{9i). Note that, this is equiva- 
lent to replace the couphngs by — > cos(6'i) cos{6j). 
The advantage of this mapping is that we can use the 
standard SC-LRPA to calculate the Curie temperature 
Tc, although the ground state is canted. Note also that, 
this mapping is exact in the two limiting cases: (i) small 
superexchange coupling and (ii) large Jaf limit where 
the canted pairs become deconnected from the system. 
As will be discussed in the following, a comparison with 
Monte Carlo results supports this procedure in the inter- 
mediate coupling regime as well. 




FIG. 2: (Color online.) Curie temperature as a function of the 
superexchange coupling strength. The density of magnetic 
impurities is set to x — 0.03 and the coupling range A varies 
from 0.2 to 0.75. Ji — Joe denotes the nearest neighbour 
coupling in the absence of superexchange. The dashed line 
shows the instability threshold of the simple (non-canted) SC- 
LRPA calculation. 

We now come back to the simple model where the 
long ranged couplings are defined by Jij — Jqc^^ . In 
Fig. [21 we have plotted the Curie temperature as a func- 
tion of Jaf for a fixed density of magnetic impurities 



and various values of the coupling range. First, as long 
as the superexchange coupling is smaller than typically 
the ferromagnetic coupling between nearest neighbors, 
(Ji = Joe~75A)^ Tc is almost insensitive to Jaf- This 
is in agreement with previous observation that in the 
diluted regime, the Curie temperature is controlled by 
couplings corresponding roughly to the average distance 
between the magnetic impurities, x~'^/^. When Jaf is in- 
creased, there is a critical value above which the ground 
state gets canted (the critical value increases with A). 
When this happens, we observe a reduction of Tc, and 
then a saturation to a finite value for strong Jaf- Let 
us discuss the limit of strong Jaf- The saturation of 
Tc corresponds to the regime where nearest neighbour 
spins are orthogonal to the other spins. The saturated 
value can be viewed as that of a system of Xeff spins 
(interacting with ferromagnetic interactions) in which all 
pairs have been removed. For example, on the fee lat- 
tice, for X — 0.05 the concentration of spins involved in 
pairs is 0.0225, or Xeff — 0.55 x. The new characteris- 
tic distance between remaining impurities is x^^j^ and 
has to be compared to the coupling range A, which con- 
trols whether the impurities "percolate" . When reducing 
A, the ratio of the saturating value over Tc{0) gets quite 
large (see Fig. [2]) because we approach the regime of non- 
"percolation" , where the remaining impurities get more 
weakly coupled. To conclude this paragraph, we observe 
a smooth crossover in the Curie temperature between 
a weak-coupling regime where Jaf has no effect and a 
strong-coupling regime where it is equivalent to remove 
all spins which are coupled by Jaf- This is true for a 
concentration of impurities small enough and we now in- 
vestigate the effect of varying the concentration. 




FIG. 3: (Color online.) Curie temperature (in units of Ji) 
as a function of the magnetic impurity concentration x for 
different values of the superexchange strength. The range of 
the couplings. A, is fixed to 0.50. 

In Fig. [31 we have plotted the Curie temperature as 
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a function of the impurity concentration, x, for a fixed 
value of the parameter A and different values of Jaf- 
In the absence of superexchange coupling, we observe a 
strong increase of Tc with x. For Jaf — 4Ji, the Curie 
temperature is strongly reduced and exhibits a maximum 
at a; « 0.15. The maximum reduces to x w 0.08 and is 
more pronounced for Jaf = 10 Ji. Above this value the 
Curie temperature decreases strongly to eventually van- 
ish at X w 0.30. The presence of a maximum can be un- 
derstood as resulting from the competition between two 
effects. As we increase the density of magnetic impuri- 
ties, X, the local fields increase (the impurities interact 
more strongly because they get closer). At the same time 
the probability of nearest neighbor spins increases also. 
Since nearest neighbor spins get canted the local fields 
they create on the other spins is reduced (and eventually 
vanish in the limit of infinite Jaf), so that the number 
of magnetically-active spins (as far as ferromagnetism is 
concerned) is effectively reduced. Note that for strong 
Jaf, one expects a site percolation which for nearest 
neighbour coupling on fee lattice occurs at xc ~ 0.20. 
Beyond this critical value, the phase should be of Neel 
type. One would need to include all fluctuations in or- 
der to calculate its critical Neel temperature, a problem 
beyond the scope of the present paper, since we focus on 
the ferromagnetic phase only. 
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FIG. 4: (Color online.) Calculated phase diagram 
(temperature-antisite concentration) for Gai-aj-^ AsyMn^^As 
indicating the predicted canted phase. The antisite concentra- 
tion y allows to tune the carrier concentration. The density of 
magnetic impurity x is fixed to 0.05. The transition temper- 
ature is calculated with modified LRPA, and compared with 
Monte Carlo (from ref. 23]), using the same set of ab-tnitto 
couplings. m(0) is the total magnetization per spin at zero 
temperature. 



III. REALISTIC COUPLINGS: THE CASE OF 
GaMnAs 

Let us now discuss the case of the widely studied di- 
luted III-V magnetic semiconductor Gai_a;Mna;As. First, 
we remind that the substitution of Ga'^"'" by Mn^+ intro- 
duces a localized spin S=5/2 and a hole in the valence- 
band (more precisely in the impurity band). 

During the molecular beam epitaxy (MBE) growth 
of the samples, there are additional defects which ap- 
pear, namely As anti-sites (Asq^^) which substitute the 
Ga sites. Arsenic antisites formation is one of the main 
mechanism for compensation in diluted magnetic semi- 
conductors. They lead to the reduction of the density of 
carriers which in turn reduces the strength of the mag- 
netic couplings and eventually the Curie temperature. In 
term of holes, Asqq^ is a double acceptor (double donor 
of electrons). If y denotes the density of As antisites 
and X the density of Mn^+ then the density of holes is 
Uh = X — 2y. The reduction of the carrier density via 
Asqj^ not only affects the long range couplings but also 
allows to "tune" the superexchange coupling22i^, and 
thus provides a way to test the ideas developped above. 
Indeed, beyond a certain concentration of Asqq^ the su- 
perexchange mechanism dominates the nearest neighbour 
coupling (see Figure 5 in [23]). This is clear from ab 
initio studies, where the couplings have been calculated 
without any adjustable parameters (Tight-Binding Lin- 
ear Muffin Tin Orbitals), and subsequently used in sev- 
eral publications^i^ii^ 



In Fig.m using SC-LRPA in the same way as described 
in section II, and using ab-initio couplings, we calcu- 
late the predicted phase diagram for Gai_a;_j,AsyMnj:As 
(temperature - density of As antisites). By energy min- 
imization, we have found a wide region of the phase di- 
agram (0.01 < y < 0.0175) where the ground state is 
not fully ferromagnetic but canted. This is the conse- 
quence of the superexchange coupling Jaf that increases 
when the antisite concentration gets larger. As far as the 
transition temperature is concerned, we stress that the 
values obtained are in very good agreement with that of 
Monte Carlo simulationsi^ Although the nature of the 
phase was not discussed in their study, this validates the 
simplified treatment of the canting of the ground-state. 
Note that, both Monte Carlo simulations and SC-LRPA 
were performed (i) with the same exchange couplings and 
(ii) with the same number of shells (approximately 20). 
Note also that the values of the Curie temperature in 
the fully polarized ferromagnetic phase {y < 0.01) are a 
little smaller here than those published in ref. [23| be- 
cause the number of shells is smaller. More importantly, 
in this study, the Curie temperature vanished abruptly 
beyond y — 0.01 (see Fig. 1 of ref. ,27]). This insta- 
bility in fact signals the occurence of a new phase, that 
we identify as a canted phase. Thus, in contrast to what 
was claimed before, the ferromagnetism survives down to 
much smaller concentrations of carriers than anticipated 
(up to y ~ 0.0175), but is non-saturated. Beyond this 
concentration, the canted phase disappears because the 
long-range couplings also become antiferromagnetic and 
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FIG. 5: (Color online.) Distribution P{9) of the canting angle 
of the nearest neighbour pairs in Gai_a;_j; AsyMn^ As for two 
different values of the concentration of As anti-sites y. The 
density of Mn^+ is s = 0.05. 



thus introduce real frustration into the system. Because 
of that, for y > 0.0175 the ground-state is expected to 
be a spin-glass. 

To characterize the ground state, we have also calcu- 
lated the total magnetization per spin m{Q) sd, T = K . 
For y < 0.01 the magnetization is to(0) = 1 (by defini- 
tion). As we enter the canted ferromagnetic phase the 
total magnetization starts to reduce significantly. For 
y = 0.0125 the magnetization is already reduced by 10% 
and for y = 0.015 it is only to(0) « 0.60 which is very 
close to the lower bound obtained by removing all pairs, 
m{0) > Xeff/x = 0.55 for x — 5%. In order to get 
an idea of how the magnetization changes from site to 
site, we have plotted the distribution of canting angles 
in Fig. [5] Note that the distribution is given for spins 
having at least one nearest neighbor. We observe for 
y — 0.0125 that approximately 40% of the spins are still 
not canted (delta peak at 6* = 0).The distribution of an- 
gles is very broad with a maximum at about 50°. On the 
other hand, for y — 0.015 we observe a strong change in 
the distribution. In the latter case, all spins are canted, 
and the distribution peaked at about 9 « 80° is more 
narrow than that of ?/ = 0.0125. 

Let us now discuss the relation between our calcula- 
tions and experimental data. It is often seen in the 
literature that the measurement of the bulk magne- 
tization (by SQUID) is different from the magnetiza- 
tion expected from the determination of the Mn den- 
sity from XRD measurement, assuming a fully polarized 
ferromagnet .^^'^^i^°'^^i^^ The direct measurement often 
leads to much smaller values. Furthermore, it is also seen 
that the magnetization strongly changes after annealing 
of as grown samples. During annealing, the magnetic im- 



purities are redistributed in the sample, which becomes 
more homogeneous.'^^ In a recent study,^^ it was shown 
that one could explain the effect of different annealing 
treatments'^'^ by the existence and the rearrangement of 
interstitial Mn defects (Mn/). Indeed Mn/ is a defect 
that preferentially sits near a Mn ion (which substitutes 
Ga) and is coupled antiferromagnetically to it . — > — This 
leads to the formation of a local singlet state for the 
dimer of Mn and therefore reduces the number of mag- 
netically active Mn, and hence the total magnetization. 
However, it is now possible to control the density of car- 
riers by chemical hydrogenation of the samplesJ^i^^ In 
this process, it is believed that the density of Mn and 
the density of defects Mn/ do not change. Therefore, if 
interstitial Mn/ defects were the main source of reduc- 
ing the bulk magnetization m(0), one would expect m(0) 
to remain the same for all these hydrogenated samples. 
This is in contradication with the measurements that in- 
dicate that the samples with lowest carrier density (in- 
sulating or very dirty metallic behavior) have a much 
smaller m(0). For instance, in Fig. 3 of ref. [s^ at small 
fields {H = 500 Oe), we observe that the magnetization 
is about two times smaller for hydrogenated sample with 
the lowest Tc, compared with the reference sample with 
no hydrogenation. This is very hard to reconcile with 
the presence of Mn/ because this would require a large 
number of such defects. Our study points out a differ- 
ent compensation mechanism, that must be at play once 
the coupling between nearest neighbours is antiferromag- 
netic (as evidenced from ab initio studies). As we said 
before, the number of pairs of nearest neighbors is large at 
X = 5%, so the reduction of m(0) is already large without 
having to invoke a large number of Mn/. We therefore 
argue that the reduction of the total magnetization is due 
to the canting of the pairs which occurs when the density 
of carriers becomes small enough, and we suggest to re- 
analyze the experimental data on the basis of the present 
work. 



IV. MODEL WITH RKKY OSCILLATIONS 

In the last paragraph we discuss another interesting 
case, the interplay between frustration resulting from the 
RKKY couplings and the superexchange. In order to 
compare with previous studies, we assume the long range 
couplings to be given by = Joe~^F(fc/?r), where 
F{kpr) — {kpr)cos{kFr)/{r/ro)^ and r = jr^ — rj\, tq 
is the nearest neighbour distance . The parameter is an 
effective Fermi vector kp which is determined by the den- 
sity of carriers. Note that this study is also motivated by 
the fact that RKKY couplings are often used to study the 
ferromagnetism in diluted magnetic semiconductors,—!^ 
although it was shown that they are inappropriatci^ 
In previous studies, it was argued, in particular, that 
the stability region for ferromagnetism was very narrow, 
upon increasing Jj^p ^ a point that was reaffirmed later 
using Monte Carlo simulations.?- We now argue that the 
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FIG. 6: (Color online.) Curie temperature for the "RKKY" 
model as a function of kp for different values of the su- 
perexchange coupling Jaf- The density of impurity is set 
to a; = 0.05 and the parameter A — 0.50. The calculations are 
performed on the fee lattice. 



stability region is in fact wider thanks to unsaturated 
phases. 

In Fig. [B] we have plotted the Curie temperature as a 
function oi kp for the RKKY-like model defined above. 
In the absence of superexchange coupling we observe that 
the Curie temperature (solid circles) exhibits a maximum 
and vanishes above kp ~ 0.60 because of the oscillations 
in the couplings. This is in agreement with ref. [2^ . 
However, when we switch on the superexchange coupling, 
we observe that the Curie temperature is reduced but 
does not vanish. This is in contrast to what was pub- 
lished previously where the ferromagnetism was appar- 
ently suppressed by the superexchange coupling. The 
reason for this discrepancy is as discussed previously, the 
occurrence of a new phase. Indeed previous calculations 
were done without including the canting. Again, the dis- 
appearance of Tc (Fig. 5 of ref. [l^) reflects only the 
change in the nature of the ground state. As we increase 
Jaf further, we observe a saturation in the Curie tem- 
perature and almost an unchanged region of stability of 
the ferromagnetic phase. These results are also in con- 
tradiction with those of ref. [1^. Indeed, it is shown 
in this paper that the ferromagnetic phase is similarly 
suppressed (see Fig. 3 of ref. [20|). The drastic reduc- 
tion of the ferromagnetic phase occurs also at low carrier 
concentration, although the frustration effects due to the 
long range couplings are very small in this region. Let us 



discuss the origin of the discrepancy. In the absence of 
Jaf, Monte Carlc^'' and SC-LRPA2^ give the same re- 
sults (see Fig. 1 and 4 of ref. [H and Fig. Ic of ref. ^26]). 
The effects of frustration are therefore properly handled 
by SC-LRPA, and thus it cannot be the source of the 
disagreement. In addition, the comparison between the 
Monte Carlo simulations of [1^ with the modified SC- 
LRPA clearly shows that the effects of the superexchange 
are also properly treated. We suggest that the discrep- 
ancy comes from the existence of the canted phase that 
was missed before. It would be of interest to clarify this 
issue. 



V. CONCLUSION 

In conclusion, we have shown that the competition 
between nearest-neighbour antiferromagnetic superex- 
change coupling and long-range ferromagnetic couplings 
gives rise to a canted ferromagnetic phase in dilute mag- 
nets. We emphasize that short-range and long-range 
competing interactions play quite different roles. In the 
later case, the oscillating tail of the RKKY interaction, 
for instance, introduces frustration at long distance and 
consequently reduces the stability region of the ferro- 
magnetic phase. In the first case, however, when su- 
perexchange is added, pairs of spins get canted but the 
stability region remains weakly affected, even at strong 
coupling. The Curie temperature of the canted phase is 
reduced simply because the canting weakens the internal 
fields. More generally, in random dilute systems, com- 
peting couplings of a range much shorter than the typi- 
cal inter-impurity distance should lead only to local spin 
reorientations, and will not affect the long-range proper- 
ties. Applying these ideas to GaMnAs, we have predicted 
the existence of a canted phase in the phase diagram. We 
have calculated its critical temperature using a modified 
local RPA approach that was found to be reliable by com- 
parison with Monte Carlo simulations using the same ab 
initio couplings. The existence of this phase provides a 
simple explanation to recent experiments in diluted mag- 
netic semiconductors, where the bulk magnetization was 
found to be smaller than the saturation value; without 
having to invoke a large number of compensating defects 
in the samples. It would be of great interest to check by 
local probes whether the ground state is indeed canted 
in the range of concentrations we have predicted. 
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